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Abstract of EP0718746 

The multiplicand (X) is multiplied (10) by a factor 
(Y) taken from a trigonometric function store (14) 
and encoded (12). The argument (theta) of the 
function is entered a binary word in the store, 
which holds only values relative to a single 
quadrant. If the input value lies outside that 
quadrant, a conversion unit (11) transforms it in 
accordance with the logic states of the most 
significant bits (BP-1,BP-2). The values of the 
trigonometric function may be obtd. by 
interpolation. 
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(54) MuKiplieur de booth pour fonctions trigonometriques 



(57) Circuit pour multiplier des donn6es conform^- 
ment a un algorithme de Booth (10) dans lequel le co- 
dage (12) des comnoandes est adapts aux caract^risti- 
ques de symetrie d'une fonction trigonometrique. On 
stocke tes valeurs de la fonction relatives d un seul qua- 
drant dans une nrtemoire (1 4) et une unite de conversion 



(1 1 ) permet d'utiliser le circuit quel que solt le quadrant. 
Les valeurs de la fonction trigonometrique peuvent etre 
obtenues par interpolation. Application aux calculs de 
transformees de Fourier ou de transformees DOT. 

Application : transformees de Fourier, transform^e 
DOT 
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Description 

L'invention concerne un circuit pour multiplier des donn6es par des valeurs d'une fonction trigonomdtrique d'une 
variable d'entr6e, le circuit comprenant : 

5 

des premiers moyens pour executor une multiplication en appliquant un algorithme de Booth, 

des seconds moyens pour gdnerer des commandes appliquees aux premiers moyens, dont des commandes de 

signe, conform6ment h I'algorithme de Booth, 

des troisi^mes moyens pour delrvrer les valeurs de la fonction Irigonometrique. 

10 

Elle concerne aussi un systdme de traitement num^rique utiltsant un tel circuit et plus particuli^rement un syst^me 
pour un calcul de transform6es de Fourier ou autres. Elle conceme aussi un systfeme de transmission num6rique 
comprenant un circuit numerique de recuperation de porteuse. 

II est courant d'utiltser des circuits pour le calcul de transformees de Fourier dans le traitement numerique de 
IS signal. Une transformde de Fourier directe et une transform^e de Fourier reciproque concement pour la premiere une 
transformation de signaux representes dans le domaine tempore I en signaux representes dans le domaine frequenciel 
et pour la seconde une transformation r6ciproque. Ces transformations qui mettent en oeuvre des fonctions exponen- 
tielles complexes impliquent la manipulation de fonctions trigonometriques. Ces manipulations necessitent de realiser 
un nombre Important de multiplications par des fonctions sinus/cosinus. 
20 It est done souhaitable d'utiliser des techniques rapides de multiplication pour diminuer les temps de calcul en y 

associant, si possible, I'utillsation d'un materiel rdduit pour timiter les coOts de mise en oeuvre. 

Une technique connue rdpondant k ces specifications est celle de I'algorithme de Booth qui consiste a grouper 
les bits des op6randes d'une multiplication par groupes de 3 bits et a generer un nombre reduit de signaux de com- 
mandes, par exemple trols commandes seulement, par groupe de trois bits : 

25 

ZERO : ne rien faire 

SGN commande le signe de I'operation ^ effectuer 

X2 : d6termine si le muttiplicande doit etre ajout6/retranch6 (selon SGN) une fols ou 2 fois aux sommes partlelles. 

30 Un muttiplieur qui met en oeuvre I'algorithme de Booth, ^ I'aide de ces trols commandes, realise des multiplications 

rapides. 

N6anmoins, des techniques permettant d'am6liorer encore ces performances sont toujours recherch6es. 

On connait le document intitule : "Parallel architecture modified Booth multiplier" de A.R. COOPER, lEE Procee- 
dings Vol. 1 35, Pt. G., N** 3, June 1988, qui d6crit un multiplieur de Booth am6lior6. 
35 En regroupant les operations select ion nees par les sorties des encodeurs de Booth et en les faisant operer en 

parall^le, il est possible de simplifier le d6roulement des operations r6alis6es selon I'algorithme de Booth. Un tel mul- 
tiplieur peut etre utilise dans tout calcul et ne tient pas compte de la particularite de I'application done des donn6es qui 
sont k multiplier entre-elles. 

Un objet de invention est de simplifier encore la construction d'un multiplieur de Booth en tenant compte de la 
40 sp6cificit6 des donn6es notamment dans I'exploitation des donn6es trigonometriques. 

Ce but est atteint avec un circuit pour multiplier des donndes par des valeurs d'une fonction trigonometrique d'une 
variable d'entr6e, tel qu'll est d6crit dans le pr6ambule, dans lequel les troisl6mes moyens dSlivrent des valeurs de la 
fonction trigonometrique sur un domaine llmlte de la variable d'entree, et en ce que le circuit comporte des moyens 
pour modifier les commandes de signe d'apr^s des etats logiques de bits de poids fort de la variable d'entree et des 
-^5 moyens de conversion de la variable d'entree pour faire que le circuit fonctionne avec toute valeur de la variable 
d'entree. 

Ainsi avantageusement, en stockant uniquement les valeurs pour un domaine Ilmlt6 de la variable d'entr6e (par 
exemple -nJ2 <e <0 pour une fonction cos 6), on reduit tres serleusement la taille de la memolre. De plus, en faisant 
que la commande de signe sort modifi6e apr6s les moyens qui g6n6rent classiquement les commandes, on realise 

so simplement avec peu de nnateriel la commande des inversions de signe qu'il est necessaire de faire pour realiser la 
multiplication mdme lorsque la valeur 9 d'entrde n'est pas dans le quadrant trigonometrique correspondant aux valeurs 
stockees. En particulier, ceci evite d'avoir k ajouter un multiplieur par - 1 pour effectuer des changements de quadrants 
ce qui consomme k la fois du temps de calcul et de la surface de circuit integre. 

Preferentiellement les troisiemes moyens dellvrent des valeurs negatives exprimdes en complement a 2. AInsi 

ss avantageusement, pour un format binaire donne (nombre de bits et repartition des poids binaires), il est possible de 
reprdsenter exactement la valeur - 1 , ce qui ne seralt pas possible avec la valeur + 1 , exprim6e en complement k deux, 
sans introduire une imprecision de calcul de 1 bit du poids le plus falble. Ceci constitue un avantage important dans 
le cas de multiplications effectuees k I'aide de I'algorithme de Booth decrit. 
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Seton un mode particulier de realisation, les moyens qui delivrent !es valours de la tonction trigonomdtrique cal- 
culent ces valeurs par interpolation autour d'un nombre Iimit6 de valeurs de fonction pr6alablement nn6moris6es. Ceci 
rdduit 6galement s^rieusement la taille de la memoire. 

Ces dlff6rents aspects de {'invention et d'autres encore seront apparents et 6Iucid6s k partir des modes de r6ali- 
5 sat ion decrlts ci-apres. 

L'inventlon sera mieux comprise ^ I'aide des figures suivantes donn^es h titre d'exemples non limitatifs qui 
repr6sentent : 

Figure 1 : un schema d'un circuit multiplieur selon invention. 

Figure 2 : un schema d'un circuit multiplieur selon invention muni de moyens pour interpoler des valeurs de fonc- 
fo tion entre des valeurs stock6es dans un premier bloc-m6moire. 

Figure 3 : un graphique d'une partie agrandie d'une fonction F(G) - cos (0). 
Figure 4 : un schema d'un circuit d'interpolation. 

La figure 1 represente un cicult multiplieur selon invention comprenant des moyens 10 MULT pour multiplier un 
multiplicande X par un multiplicateur Y conformement a I'algorithme de Booth, des moyens 12 ENCOD pour g^ndrer 
IS des commandos conformement a cet algorithme de Booth afin de permettre aux moyens 10 d'operer, et des moyens 
14 MEM pour d6livrer les valeurs du multiplicateur Y au multiplieur 10. 

L'algorithme de Booth consiste ^ ajouter un z6ro en bit de potds faible puis ^ regrouper les bits du multiplicateur, 
exprimS en complement k 2, en groupes de trois bits consecutifs ^ partir du bit de poids faible ajoute et h en deduire 
tea operations ^ effectuer. Deux groupes cons6cutrfs se chevauchent par un bit. Les trois bits d'un groupe Imposent 
^0 le type d'ope rations k effectuer : 



Code 


Operation 


(001; 010) 
(000; 111) 
(Oil) 

(101, 110) 
(100) 


- additionner le multiplicande k la somme partielle 

- ne rien faire 

- additionner deux fois le multiplicande k la somme partielle 

- soustraire le multiplicande de la somme partielle 

- soustraire 2 fois le multiplicande de la somme partielle. 



Ces cinq commandos peuvent etre ramenees a trois commandos soulemont : 

- ZERO - ne rien faire 

- SGN - commando le eigne de I'op^ration k effectuer 

- X 2 - determine si le multiplicande doit etre ajoute/rotranche (selon SG N) une fois ou 2 fois de la somme partielle. 

II y a done autant de fois les trois commandos prececontes qu'il y a do triplets de bits dans le multiplicateur Y Si 
celui-ct contient b bits, le nombre de triplets est egal k rentier lo plus proche soit de b/2 soit de b/2 + 1 selon la parit6 de b. 

L'invontion concerne des modifications de ta commande du signe SGN delivree normalement par un encodeur de 
Booth (12). 

Lo multiplieur 10 de Booth revolt lo multiplicande X et plusieurs lots des trois commandos pr6c6dontos. Le multi- 
plicateur Y est cod6 par les moyens de codage 12 ENCOD pour d6livrer la pluralit6 do commandos ZERO, SGN et 
X2. Dans l'invontion, le multiplicateur Y represente des valeurs d'une fonction trigonomdtrique, par exomple cos 9. II 
est obtenu en adressant des moyens de memorisation 14. par exemple une m6moire MEM, qui ont pr6atablement 6t6 
charges par des valeurs de la fonction trigonometriquo. L'argumont G do la fonction trigonometrique est fourni en entree 
du circuit selon un mot binaire Bp.i-.-.BQ. En utilisant les propri6t6s de sym6trie de la fonction cos 0, il est possible de 
ne stocker dans la memoire 14 que les valeurs cos 9 corrospondant k 7i/2 <9 <n, c'est-^-dire des valeurs negatives de 
cos 0. Les deux bits do poids fort Bp.i ot Bp.2 permettent do deduire les valeurs de cos 0 pour los autres valeurs de 9. 

Selon l'invontion, on controto la commande de signe SGN de I'encodeur de Booth. Une porte 18 re9oit Bp.^ et Bp.2 
et delivro un signal CD. Une famille 19 de portes OU - exclusif regoit le signal CD et les signaux SGN issus de chaque 
triplet de bits de Y En sortie, on obtient de nouveaux signaux do signe SGN' qui sent utilises pour s6lectionner le type 
d'oporation k effectuer dans lo multiplieur de Booth 10. Le mot binaire Bp.^.. Bq est cod6 en complement k 2. Chaque 
nouvelle commande de signe SGN' est alors deduite de la table I de verite suivante : 



TABLE I 





Bp.2 


CD 


SGN 


SGN' 


0 


1 


0 


0 


0 


0 


0 


1 


0 


1 



3 



EP 0 718 746 A1 



TABLE I (suite) 



1 


0 


0 


1 


1 


1 


1 


1 


1 


0 



5 

Preferentiellement, la m6moire 14 stocke des valeurs negatives de cos 9 done avec jt/2 <0 <7t ou n<Q <3n/2. 

Lorsque la valeur 6 qui arrive en entree du circuit n'est pas dans le quadrant correspondant aux valeurs stock^es 
dans la memoire, alors une unrtd de conversion 11 transforme la valeur 0 d'entr^e d'apr^s les etats logiques des bits 
de poids 6lev6 Bp.^, Bp.2. Selon le contenu de !a rT>6moire 14 (sinus ou cosinus), la conversion va consister k convertir 
10 6 en n/2 - 9 et a modifier eventuellement le signe selon le quadrant concern^. 

Pour r^uire encore les moyens mat6riels, les moyens 14 qui ddlivrent la fonction trigononn6trique peuvent op6rer 
par interpolation autour d'un nombre limits de valeurs de la fonction trigonom^trique prealablement nnemorisees (figure 
2). 

Le mot binaire Bp.^... Bq est alors separ6 en trots champs : 

15 

un premier champ I form^ par des bits de poids faible du mot binaire, 
un second champ II fomrid par des bits medians du mot binaire, 

un troisidme champ III, form6 par les bits Bp.-,... Bp.g, qui sert k contrdler les commandes de signe comme expos6 
auparavant. 

20 

Les premier et second champs I et II servent k adresser (figure 2) une unite de conversion 1 1 CON V qui transforme 
les valeurs d'entree 9 en valeur de 9 situ6e dans le quadrant pour lequel on a memorise prealablement un nombre 
limits de valeurs de la fonction trigonom^trique dans un premier bloc-m6moire 142. Les valeurs stock6es peuvent ^tre 
prealablement corrigees par des coefficients correctlfs spdcifiques destines k reduire des erreurs engendrees par la 

25 m6thode d'interpolation utilis6e. 

L'unite de conversbn adresse aussi un second bloc-memoire 144 dans lequel on a memorise prealablement des 
coefficients d'interpolation calcu!6s pour les valeurs de la fonction trigonom6trique qui ont 6t6 stock6es dans le premier 
bloc-memoire. Les coefficients d'interpolation peuvent etre une derivee vraie ou une deriv6e approximative formee par 
exemple par ia pente moyenne du segment joignant deux valeurs stock6es cons6cutives de la fonction. 

30 Ainsi en adressant les deux blocs-memoires avec le champ median II (avec ou sans conversion selon le cas), par 

exemple on peut lire la valeur de la fonction cos 9 (avec corrections) et celle des coefficients d'interpolation pour cette 
valeur 9. 

Le champ I de bits de poids falbles est alors exploit6 pour calculer une quantity interpol6e li^e k I'incr^ment fourni 
par le champ I. 

3S La figure 3 repr^sente une partie d'une fonction F(9) = cos(9) connue par un nombre limits de valeurs F(9nk) = cos 

(^nk) [point A], F(9(„^.ijk) = cos(9(^ij;J [point B], pour des valeurs e^k, 9(n+i)k variable 9. Sur le graphique n est 
rindice courant des points A, B..., avec 0 ^n<^-^ , k 6tant un nombre entier d6terminant le facteur d'interpolation, avec 

Lorsqu'une valeur quelconque G^j^^j de la variable est pr6sent6e en entrde du circuit, celui-ci calcule, par interpo- 
40 lation, une valeur interpolee proche de la valeur f{Qnf^Ji [point C]. Sous cette forme, I'indice i est un indice courant 
variant de 0 ^ k-1 afin de diff6rencier k valeurs intermedial res entre G^k et 9(^^.1)1^. 

Pour repr^senter N valeurs de la variable 0, il faut log2N bits. Une valeur quelconque de la variable est 
representee par un mot binaire pour lequel on ignore les deux bits de poids forts du debut pour ne consid6rer que les 
valeurs appartenant a un seu I quadrant. Le mot binaire 1 3 (figure 2) comprend un premier champ 1 3A comportant des 
4S bits de poids fort et un second champ 1 3B comportant des bits de poids faible. 

Selon rinventton, 1' Interpolation est effectuee en exploitant un nombre N/k de valeurs connues de la fonction. Ces 
N/k valeurs connues sont preferentiellement stockees prealablement dans un bloc-memoire 142. Les valeurs connues 
F(9p,)() sont celles qui correspondent aux log2N/k bits de poids fort de la variable d'entree. 

Pour les N/k valeurs connues de la fonction (points A, B, etc), selon invention, on stocke egalement. preferen- 
50 tiellement dans des moyens de memorisation 144, N/k coefficients d'interpolation C(9„j()F Ceux-cl peuvent etre des 
valeurs de derivee F'(9„^) de la fonction F(9n,^). Ces valeurs de derivee peuvent etre des valeurs moyennes de derivee 
calculees autour de chaque point connu A, B, etc. Preferentiellement, ces valeurs de coefficient d'interpolation sont 
egales k la pente des segments reliant deux points connus consecutifs. 

Ces valeurs de coefficient d'interpolation peuvent egalement dtre formees par des valeurs vrales de derivee cal- 
55 cuiees aux points connus A, B. etc. 

Pour limiter les Imprecislons, on stocke dans le bloc-memoire 142 des valeurs connues de la fonction prealable- 
ment corrigees par des termes correctlfs. Le bloc-memoire stocke done des valeurs f(%i) + cor(Qf^y). Le terme cor 

est prealablement calcule pour chaque point connu pour minimiser I'erreur quadratique moyenne entre les valeurs 
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interpolees et la fonction Irtgonomdtrique parfaite. 

Lorsqu'une valeur d'entr^e est prSsentSe d I'entr^e du circuit, le champ 13A de bits de poids fort adresse 
(connexion 15) le btoc-memoire 142 et les moyens de memorisation 144 : 

5 - le bloc-memoire delivre F(e„,^) + cor(9„^) constituant la valeur approchee de la valeur interpolee k determiner 
(connexion 17), 

les moyens de memorisation delivrent les coefficients d'interpotation C'{%^^) (connexion 1 6). 

Conslderons le cas ou la fonction F(6nk) fonction cosinus et ou les valeurs de coefficients d'interpolation 

10 sont d6duites de la pente d'un segment y a x+b. Ce segment passe au voisinage des deux points connus A et B 
pour indiquer que des coefficients de corrections interviennent dans le calcul des valeurs interpolees. 
Une valeur approch6e y de cos (Spk+j) [point D] est alors telle que : 

cos 0-,t^|,-cos 0 . 

y= e -e ^-'^^ ^^^^ ^-f< 

IS ^(nk+k) ^nk 

avec 

0<i<k-^ 0<n<^-1 
Aprfes des manipulations trigonom6triques connues, on obtlent : 

V - r o (2n+l)kn ^. kn] i ^ 2n . 

Y - -2 sxn -r .sin — - — + cos — - nk 

N N \ k N 

2S Pour obtenir une valeur interpol6e cos (©nk+j) plus proche de la valeur exacte cos (6nk+i)' introduit un terme de 

correction cor [%^). Cette valeur interpolee est situee entre les ordonnees des points D et C. Chaque terme de cor- 
rection cor (e^k) ®st determine en rendant minimum la somme Xfcos 6^;^+/ - cos Onk+t)^ 9^^^ 0 ^ i < k-1 . 
La valeur approch6e cos (e„^^^) s'6crit alors : 



30 



cSs (e„,.,) = j-| sin (i^./cnj.sin (^) } i ♦ cos nk) * cor <e„,) 



Le circuit comprend des moyens de calcul 145, par exemple un multiplieur/accumulateur qui revolt 0*(%i^) et le 
35 champ 1 3B forme des bits de poids faible de la valeur d'entree 9nk+i- 

Dans une premiere etape, les moyens de calcul 145 calculent le resultat partiel : 

et dans une ^econde 6tape, ils ajoutent ce r6suttat partiel k la valeur approch6e f(Qn\) + ^o^{%k) P^^^ d^livrer la valeur 
40 interpolee : fi^nf^+H soit : 

Ainsi, il est possible de r6duire la tailte de la seconde m6moire en ne stockant qu'un nombre limits de valeurs sur 
un nombre r6duit de bits tout en obtenant une precision 6lev6e sur le r6sultat, ce qui est particuli^rement int6ressant 
45 dans la realisation de circuits integres. 

Dans le cas ou ils delivrent des valeurs interpolees, les troisidmes moyens comprennent : 

le premier bloc-memoire 142 pour m6moriser un nombre iimite de valeurs negatives de la fonction correspondant 
au second champ de bits, a des adresses issues du second champ de bits, 
50 - le second bloc-memoire 144 pour memoriser des coefficients d'interpolation associes aux valeurs memorisees 
dans le premier bloc-memoire, a des adresses issues du second champ de bits, 
des moyens d'interpolation 

pour lire, dans le premier bloc-memoire, une valeur approchee de la fonction en I'adressant par le second 
55 champ de bits de la variable d'entr6e, 

pour lire, dans le second bloc-memoire, le coefficient d'interpolation en I'adressant par le second champ de 
bits de la variable d'entree, 

pour calculer une quantite interpolee en multipliant le premier champ de bits de la variable d'entree par le 
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coefficient d'interpolation adressd, la quantite interpol6e 6tant totalisee avec la valeur approchde pour donner 
la valeur interpol6e. 

Le circuit peut dtre utilisd dans un systdme de traitement num^rique de donn^es pour eftectuer des catculs de 
5 transformdes de Fourier, de transformees en cosinus discretes, ou d'autres transformees faisant appel ^ des fonctions 
trigonom6trlques. 

Le circuit peut ausst etre utilise dans un syst^me de transmission numdrique dans lequel le circuit selon invention 
entre dans la realisation d'un circuit num6rique de r6cup6ration de porteuse. Dans ce cas, la variable d'entr6e est une 
erreur de phase 0 k partir de laquelle le circuit multiplie des valeurs cos 6 et sin 6 par un multipltcande forme par un 
10 symbole qui doit subir une correction de phase pour synchroniser la r^cup^ratbn de porteuse. 



Revendications 

'5 1. Circuit pour multiplier des donn§es par des valeurs d'unefonctiontrigonometriqued*une variable d'entree, le circuit 
comprenant : 

des premiers moyens pour executer une multiplication en appliquant un algorithme de Booth, 
des seconds moyens pour g6n6rer des commandes appliqu6es aux premiers moyens, dont des commandes 
20 de signe, conformement k I'algorithme de Booth, 

des troisiemes moyens pour delivrer les valeurs de la fonction trigonometrique, 

caracterise en ce que les troisiemes moyens delivrent des valeurs de la fonction trigonometrique sur un domaine 
limite de la variable d'entree, et en ce que le circuit comporte des moyens pour modifier les commandes de signe 
^5 d'apres des etats logiques de bits de poids fort de la variable d'entree et des moyens de conversion de la variable 

d'entr6e pour faire que le circuit fonctionne avec toute valeur de la variable d'entr6e. 

2. Circuit selon la revendication 1 caract6ris6 en ce que les troisiemes moyens d6livrent des valeurs negatives de 
la fonction trigonometrique. 

30 

3. Circuit selon les revendications 1 ou 2 caract6ris6 en ce que les commandes de signe sont inversees lorsque les 
deux bits de poids fort de la variable d'entree sont identiques. 

4. Circuit selon la revendication 3 caract6ris6 en ce que les troisifemes moyens calculent la valeur de la fonction 
35 trigonometrique par interpolation autour d'un nombre llmite de valeurs de fonction pr^alablement memor isees. 

5. Circuit selon la revendication 4 caracterise en ce que la variable d'entree etant representee par un premier champ 
de bits de poids faible, par un second champ de bits de poids interm6dlaire, et par un troisi6me champ de bits de 
poids fort, les troisiemes moyens calculent, par interpolation, une valeur interpolee de la fonction trigonometrique 

40 pour la valeur de la variable d'entr6e, les troisiemes moyens comprenant : 

un premier bloc-m6moire pour m6moriser un nombre limite de valeurs negatives de la fonction correspondant 
au second champ de bits, ^ des adresses issues du second champ de bits, 

un second bloc-memoire pour memoriser des coefficients d'interpolation associes aux valeurs memorisees 
4S dans le premier bloc-memoire, a des adresses issues du second champ de bits, 

des moyens d'interpolation 

pour lire, dans le premier bloc-memoire, une valeur approchee de la fonction en I'adressant par le second 
champ de bits de la variable d'entree, 
50 - pour lire, dans !e second bloc-memoire, le coefficient d'interpolation en I'adressant par le second champ 

de bits de la variable d'entree, 

pour calculer une quantite interpolee en muttipliant le premier champ de bits de la variable d'entree par 
le coefficient d'interpolation adresse, la quantite interpolee etant totalisee avec la valeur approchee pour 
donner la valeur interpolee. 

55 

6. Systeme de traitement num6rique de donn6es caracterise en ce qu'il comprend au moins un circuit selon une des 
revendications 1 a 5 pour un calcul de transformee de Fourier ou de transformee en cosinus discrete ou de trans- 
form6e faisant appel ^ une fonction trigonometrique. 
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7. Sysleme de transmission numerique caracterise en ce qu'il comprend au moins un circuit selon une des revendi- 
cations 1 k 5 pour r^aliser un circuit numerique de recuperation de porteuse. 
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